Confined space and effective interactions of multiple self-avoiding chains 
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We study the link between three seeming-disparate cases of self-avoiding polymers: strongly over- 
lapping multiple chains in dilute solution, chains under spherical confinement, and the onset of 
semi-dilute solutions. Our main result is that the free energy for overlapping n chains is indepen- 
dent of chain length and scales as n^''*, slowly crossing over to , as n increases. For strongly 
confined polymers inside a spherical cavity, we show that rearranging the chains does not cost an 
additional free energy. Our results imply that, during cell cycle, global reorganization of eukaryotic 
chromosomes in a large cell nucleus could be readily achieved. 



How are polymers organized inside a confined space? 
What is the free energy barrier to the overlapping of two 
or more chains in the absence or presence of confine- 
ment? Is the barrier higher for longer chains? Recently, 
there has been renewed interest in the problem of con- 
fined polymers P, |3j S 01' because of its relevance to 
such biological processes as DNA packaging in a virus Q 
and organization and segregation of chromosomes in bac- 
teria [6[. In eukaryotes, multiple chromosomes are en- 
capsulated in a cell nucleus, i.e., in a dimension many 
times smaller than their natural sizes. Although a single 
cell can contain as many as billions of basepairs of DNA 
(e.g., human), the dimension > 10 /im of cell nuclei im- 
plies that the volume fraction of total amount of DNA 
is typically much less than that of virus or bacteria, and 
the eukaryotic chromosomes in a spherical volume may 
be considered as a semi-dilute polymer solution [l[. Here, 
one of the key issues is spatial organization of chromo- 
somes. The emerging view is that they are compart- 
mentalized and occupy discrete "territories" inside a nu- 
cleus 0. However, chromosomes should also be able to 
mix, when necessary (e.g., recombination), and, indeed, 
they do Q. Unfortunately, despite its importance, lit- 
tle is known about how self-avoiding polymers interact 
and are organized in a confined space, whereas the ef- 
fective interaction of polymers in dilute (bulk) solution 
has been well studied |9|. Although these two subjects, 
namely, the effective interaction in dilute solution and 
polymers under confinement, have been regarded a dis- 
parate grouping so far, as we shall show below, there is 
a close connection between the two. 

The first attempt to characterize the effective inter- 
action between two polymeric coils in dilute solution, 
each carrying N monomers, was made as early as in 
1950 by Flory and Krigbaum [lo| . The main conclu- 
sion of their "mean-field" approach is that the overlap- 
ping free energy, i.e., the free energy cost for bringing 
two chains in a volume explored by each chain, scales as 
I3Tfk iVs > 1 (Throughout this paper, (3 = l/ksT, 
where ks is the Boltzmanir constant and T the abso- 
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FIG. 1: (Color online) (a) Strongly overlapping multiple 
chains with excluded volume in dilute solution, (b) Onset 
of semi-dilute regime as "stacking spheres." Above the on- 
set concentration, each sphere is independent of one another. 
This allows one to map the overlapping chains in (a) onto 
those confined in a single sphere in (b) under equivalent con- 
ditions - the same chain (sphere) size and the same monomer 
concentration - and thus to estimate the overlapping free en- 
ergy in (a). 



lute temperature.) The repulsion between long chains 
is so strong that they should behave as mutually im- 
penetrable hard spheres. Partly due to its simplicity, 
this picture satisfied scientists for about three decades. 
Then, the more careful scaling analysis of Grosberg et 
al. [llj , which takes into account monomer-density corre- 



lations [1^, [13| , showed that the overlapping free energy 
is of the order of fc^T and is asymptotically indepen- 
dent of N. This counter-intuitive result has since been 



confirmed by renormalization group calculations 1J| and 



numerical simulations 15, 1 



The purpose of this Letter is to unravel the link be- 
tween a single chain under spherical confinement, a sys- 
tem of multiple, strongly overlapping chains in dilute 
solution, and the onset of semidilute solution [Figs. [T] 
and [2], using scaling arguments and molecular dynamics 
simulations. The polymers we consider here are flexible 
chains with excluded volume, unless otherwise stated. An 
important corollary of our analysis is a non-trivial gener- 
alization of the aforementioned two-chain result by Gros- 
berg et al. to the case of an arbitrary number of chains 
[see Fig.[lja)] in a wide range of monomer concentrations. 
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It also illuminates the weak dependence of the effective 
interaction on polymer concentration studied recently by 
Louis et al. (l6[ using numerical simulations. Biological 
implications of our results will also be briefly discussed. 

Consider a polymer solution, where its constituent 
chains start to touch one another [Fig. [IJb)]. At this 
onset of the so-called semi-dilute regime, the solution 
can be viewed as stacking of imaginary spheres, where 
the size of each sphere, containing a single chain, is ^Rg, 
I.e., the radius of gyration of each chain in the dilute 
regime. Using this "stacking sphere" picture, we can, in 
fact, apply the physics of semidilute solutions to estimate 
the interaction free energy of multiple chains in dilute so- 
lution [Fig. [Ha)]. The essence of our analysis is that, 
as the monomer concentration increases above the onset 
concentration, each sphere shown in Fig. [U^b) behaves 
as an independent replica of the neighboring ones (see 
below) . 

To this end, we first consider a long chain formed by 
N monomers, compressed inside a spherical cavity of di- 
ameter D < Rg. In the case of an ideal chain, it is well 
known that the confinement free energy can be obtained 
by a random-walk (RW) analysis, where the chain can 
be divided into "independent" subchains that start from 
one point on the wall and reach another via a random 
walk [17|. This is reasonable since, at each "collision" 
on the wall, the chain loses its memory of chain con- 
nectivity [18|. Since the length of a subchain scales as 
-^sub (D/a)'^, where a is the monomer size, the con- 
finement free energy (i.e., entropy loss of order 1 times 
—kBT) can be straightforwardly obtained by counting 
the total number of independent subchains as follows [17| 
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In the case of a strongly confined chain with ex- 
cluded volume [or a confined "self-avoiding walk (SAW)"] 
[Fig. (Ha)], the major source of increase of free energy is 
the collisions between monomers along the chain. (This 
allows us to approximate the spheres in Fig. [H^b) as in- 
dependent subsystems.) The confinement free energy 
can be obtained by mapping the chain onto an equiv- 
alent semidilute solution: In the semidilute regime, this 
monomer-monomer contact probability within the sphere 
of volume V ~ can be estimated using the des 
Cloizeaux exponent l/(3z^ — 1) The resulting free 

energy is 



(2) 



where ~ 3/5 is the Flory exponent [Ul, ll7|. Note that, 
as expected, J^saw grows with N faster than linearly, 
namely, faster than that of the corresponding ideal chain 
(Eq. [T]). Although useful to estimate the free energy. 





FIG. 2: (Color online) A self-avoiding chain under spherical 
confinement, (a) Stacking blobs in a moderately-strongly con- 
fined self-avoiding chain, i.e., -D ^ ^ 3> a. (b) Mixing and 
de-mixing of chains within a confined sphere. 



this particle picture provides little insight into the spa- 
tial organization of the confined chain. A more intuitive 
approach is that of stacking "blobs" [Fig.dja)] [i,[li,|23], 
where the size of each blob is the correlation length 
^ a{D^ /a'^ NY^^^'^^^'f in an equivalent semidilute solu- 
tion. We note that Eq.[2]is then self-consistently restored 
by the "ksT per blob" ansatz, i.e., P^saw ~ D^/£,^ 

Within the stacking-blob picture (not to be confused 
with the stacking- sphere picture), we can view a long 
chain with excluded volume confined in a spherical 
cavity as a system of overlapping multiple chains. The 
basic idea is similar to that of the RW analysis described 
above: Each "independent" subchain is a series of 
connected blobs [19|, which percolates and connects 
two points on the confining walls as illustrated by the 
red blobs in Fig. [2Ka). The number of monomers per 
independent subchain, A^subi can be obtained as follows. 

(i) Weakly confined regime {D > ^ ^ a): In polymer 
solution, this is analogous to the onset of the semidilute 
regime, where the chains start to contact one another: 
i^Rg^D, and, thus, A^sub ^ {D/af''' [Fig. [^b)] . 

(ii) Moderately-strongly confined regime (-D ^ ^ ^ 
a): This is where both the semidilute regime and the 
stacking-blob picture in Fig. [UJa) apply. In other words, 
within the correlation length ^, the chain conformation 
is that of SAW, but the global conformation of each in- 
dependent chain is described by the RW of blobs. Thus, 
iVsub ^ {(./aY'^D/af ~ (7^/a)i/(3.-i)^(2.-i)/(3.-i)_ 

(iii) Concentrated regime (D ^ ^ ~ a): In this regime, 
excluded volume is screened at all length scales beyond 
a [20]. The RW analysis then leads to A^^ub ^ (D/af. 

Using the expressions of A"sub m (i) and (ii) (the two 
regimes of our main interest), we can rewrite Eq. [2] in 
terms of the number of independent subchains, n = 
N/Nsuh, to estimate the free energy cost for over- 
lapping an arbitrary number (n) of chains. We find 
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(^ > e > a) 
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In fact, this is also the free energy cost for bringing 
n self-avoiding chains (each carrying N monomers) in 
dilute solution to an imaginary sphere of volume Vg 
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each chain would explore otherwise [Fig. [DJa)]. Here, 
Vg ^ Rg ^ {aN")^, a condition compatible with (i). 
In higher monomer concentrations compatible with (ii), 
the chain size is reduced by "screening" effects , and, 
thus, Vg ~ (a 'n}~^'^ N'^Y , which explains the larger ex- 
ponent in Eq. [H This rapid increase of with n, as 
evidenced in our blob picture, implies a strong repulsion 
between two spheres of dense chains, reminiscent of that 
between two star polymers (see Discussion for interesting 
consequences on chromosome organization). 

At first glance, this A^sub-independence is surprising - 
the increase of free energy due to chain overlapping is 
independent of the chain length, where its special case 



for n = 2 explains the results by Grosberg et al. 
This is a natural consequence of the functional form of 
Eq. [21 Furthermore, the exponent 9/4 is identical to 
the osmotic pressure exponent in the semidilute regime. 
This is not accidental: The monomer density (times a^) is 
now translated into n, which should correctly reflect the 
stronger monomer-density correlation in the semidilute 
regime [Ull^. 

Importantly, the scaling form in Eq.[2]also implies that 
the overlapping free energy is invariant under rearrange- 
ment of the chains as illustrated in Fig. [21(b), where each 
segregated chain occupies a smaller volume of linear di- 
mension d ~ {D^ /ny/^ . The total free energy after seg- 
regation is the same as before: 
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where 



Also, note that Eq s . [H and [4] ca n explain the recent sim- 
ulation results by Louis et al. [16| . Their main conclusion 
is that the two-chain interaction in a polymer solution is 
« 2kBT and independent of the monomer concentration 
c for c up to the overlap concentration c* , increasing only 
slightly with c above c*. This weak c dependence is al- 
ready apparent from Eq. [31and[31 In the scaling regime, 
the ratio of entropy loss in the dilute regime and the 
semidilute regime is only 2^/2^^'^ k 1.7, and, thus, we 
predict the two-chain interaction increases only up to 3- 
4 ksT from 2kBT. 

To augment our scaling analysis, we also performed 
Molecular Dynamics simulations using ESPResSo 21 1. 
In our simulations, the polymer chains are represented 
as a bead-spring model. We chose a chain consisting 
of A^ = 1000 beads connected by FENE (spring) bonds 
with a purely repulsive Lennard-Jones (LJ) potential 
for excluded-volume interactions (monomer- monomer as 
well as monomer- wall) . The basic length scale in the sim- 
ulation is the bead diameter a, energies are measured in 
units of the LJ interaction energy e at distance a, and 
the mass of a bead m is the mass unit. Consequently, 
time is measured in units of the LJ time t^j = a^J m/e. 
The FENE bond constant was chosen to be lOe and the 



maximal elongation of a bond as 2a. We used a standard 
velocity Verlet algorithm with a timestep of O.OItlj to 
propagate the system, and a Langevin thermostat with 
friction constant 7 = r^J to keep the system at constant 
temperature T — e/Zcs. 

We ran 40 sets of simulations to test a wide degree 
of confinement, namely, the radius of confining sphere 
R ^ D/2 hes in the range i?(,/4 < R < Rg (where 
Rg « 28.3 for A^ = 1000). We started each set of sim- 
ulations by generating a SAW in a large spherical cav- 
ity whose size was several times the Rg of the chain. 
Then, the size of the confining sphere was gradually re- 
duced, forcing the chain radially inward, until it reached 
the target value oi D (< Rg). Thereafter, we recorded 
2 X 10^ conformations and associated parameters at ev- 
ery IOtlj. To check the consistency of our simulations, 
we first computed the pressure, p, of the confined chain, 
which is related to the free energy of confinement by 
p = -^^^ or p ~ (j)TsAw/N, where F ~ DMs the 
volume of sphere and cj) ~ Na^ /D^ the volume fraction 
of the chain. Since Eq. [21 can also be expressed in terms 
of (j) as -FsAW ~ N(j)^^^^'^^^\ we obtain an equivalent 
scaling for the pressure p ~ (^3i//(3i'-i) ^ (j)^/^. This re- 
lation is indeed confirmed by our simulations (data not 
shown), which lead to p ~ ^2.27±o.02 ^^j. ^ ^ q g^^^^ 

are in excellent agreement with the recent Monte Carlo 
simulations by Cacciuto and Luijten 

Next, we tested the validity of the stacking blob hy- 
pothesis illustrated in Fig. ^a,). If the blobs can be 
considered as impenetrable hard spheres of diameter ^ 
stacked together, their radial monomer density profile in 
the spherical cavity tends to be uniform, except within a 
length of order ^ from the wall. Our simulation confirms 
this: In Fig.[3^a), we show 40 radial density profiles, p{r), 
from our simulations, which have been normalized such 
that Jf^ pir)dr — 1. Note that p{r) is indeed constant for 
a significant range of r. This feature is particularly pro- 
nounced when we compare it with that of an ideal chain 
p.^{r) = (l/27ri?)sin2(7rr/i?)/r2 [see Fig. [HJa) inset for 
the case of i? = 18.0 = O.GARg]. The decay near r = R 
manifests the existence of a depletion layer of length scale 
^ from the confining wall. One interesting observation is 
that, at high volume fraction (0 > 0.15), the density be- 
comes oscillatory near the wall. This means that the wall 
collaborates in enhancing the ordering - it is the signa- 
ture of the crossover from the semidilute regime (^ ^ a) 
to the concentrated regime (^ ~ a) mentioned above. In- 
deed, we simulated a system of A = 1000 hardspheres 
by removing the bonds between monomers in our simu- 
lations, and observed the same oscillations in p(r) (data 
not shown in Fig. [3]). 

Our final and the most important test concerns the 
view of a single confined chain with excluded volume 
as significantly overlapping multiple, independent sub- 
chains. If this view is correct, the average internal dis- 
tance between a pair of monomers i and j, R{\i — j\) = 
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FIG. 3: (Color online) Independent-subchain analysis, (a) 
Monomer densities tend to be uniform for a sizable range of 
\i — j\. The liquid-like oscillatory behavior near the wall at 
high monomer concentrations (0 > 0.4) illustrates how the 
wall cooperates in enhancing the ordering, (b) R{\i — j\) vs. 
\i — j\. The inset shows the internal distance R{\i — j\) be- 
comes saturated beyond |i — ,7| ~ A^sub, due to the "reflect- 
ing" wall. At low concentrations (weakly confined chains, i.e., 
^ ^ D, represented by blue and green curves), the data tend 
to collapse when \i — j\ is rescaled by Nguh — (R/a)^^". For 
the strongly confined case Z) S> ^ 2> a, however, the correct 
rescaling factor is (ii) N,ub = (i?/a)^/'='"-i)iV'2— i)/(3''-i) . 
These results support the independent-chain assumption. 



J I ) , should increase as their contour distance in- 
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creases up to its maximum value 
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Beyond 



i — j'l = iVsub, however, the monomers are independent, 
and their average distance is constant i?max(|i — J I) ~ R 
[see Fig. [3Ub) inset]. In other words, for any chain 
length N and a confining sphere radius i?, the reduced 
internal distance curves, R{\i — j\)/R vs \i — j\/Nsnhi 
should collapse onto each other. In Fig. [2Ib), we verify 
the independent-chain hypothesis, where we rescale the 
whole set of internal distance curves shown in the inset 
with the two different expressions of iVgub^ (i) Nsnh — 
[D/afl" and (ii) A^^ub = i?i/(3^-i)Ar(2i'-i)/(3i^-i). Note 
that, for weakly confined chains (represented in blue and 
green), the curves collapse nicely when rescaled by (i) 
[as well as by (ii) for the the size of chain N = 1000 we 
simulated] . As the degree of confinement becomes strong 



(e.g., red curves), rescaling by (ii) produces a much bet- 
ter result, supporting our argument above. In both cases, 
we note that the plateau starts at |i — j|/-^sub ~ 1- 

Our results have implications for spatial organization 
of chromosomes inside a eukaryotic cell nucleus. While 
the eukaryotic chromosomes have several levels of high- 
order structures, the length scale that characterizes the 
dsDNA and chromatin fiber is typically < 100 nm, much 
smaller than the size of the cell nucleus > 10 jim. 
It is thus conceivable that the global organization of 
these molecules in confinement will not sensitively reflect 
molecular details (e.g., structure of chromatin fiber), and, 
importantly, our analysis in Fig.[2ljb) implies that the free 
energy cost for global chromosome reorganization is low. 
On the other hand, once territories are formed (by pro- 
teins and DNA cross-linkings), the structured, compact 
and segregated chromosomes, are likely to exclude each 
other, since the overlapping of tightly packed blobs (such 
as star polymers) is highly costly (JF^) for large n, as 
mentioned earlier. Further consideration of other factors 
including the role of chain stiffness and the geometry of 
confined space is certainly warranted. 
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